Adiabatic refocusing of nuclear spins in Tm 3+ :YAG 
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We show that the optical absorbance detection of nuclear spin echo gives direct access to the spin 
concentration, unlike most coherent signal detection techniques where the signal intensity/amplitude 
is difficult to connect experimentally with the spin concentration. This way we measure the spin 
refocusing efficiency in a crystal of Tm 3+ :YAG. Given the large inhomogeneous broadening of the 
spin transition in this material, rephasing the spins with the usual hard pulse procedure would 
require excessively high radiofrequency power. Instead we resort to an adiabatic pulse sequence 
that perfectly returns the spins to their initial common orientation, at low power cost. 

PACS numbers: 42.50.Ex,42.50.Gy,42.50.Md,82.56. Jn,03.65.Wj 



I. INTRODUCTION 



II. ADIABATIC SPIN REFOCUSING 



The optical detection of magnetic resonance has been 
investigated for a long time, starting with atomic va- 
pors wj^- In solids, the earliest studies were reported in 
ruby [4j. Since the emergence of tunable lasers, special 
attention has been paid to the Raman heterodyne opti- 
cal detection of coherent transient spin processes Q . De- 
tection of coherent magnetic resonance features through 
optical fluorescence was examined even earlier [f| . In the 
latter class of processes, a tt/2 RF pulse is used to convert 
the spin coherence into a level population that can emit 
light by fuorescence. A similar approach was also used 
to detect spin echoes through optical absorbance 0- 

Optically detected magnetic resonance is enjoying re- 
newed interest as research on quantum memory for light 
becomes more popular. Most quantum storage schemes 
rely on the conversion of optical atomic coherence into 
long lifetime spin coherence and back [cj-[l^. Echo tech- 
niques are needed to compensate for the spin inhomoge- 
neous phase shift [H, [H[, i-e. to refocus the spins. 

In the present paper our aim is twofold. First we ex- 
amine the potential of adiabatic rapid passage (ARP) for 
spin refocusing, when hard pulses would require excessive 
RF power. Second, placing this research in the context of 
quantum memory, we explore the specificity of magnetic 
resonance detection through the absorption of light. We 
show that this approach may convey more information 
than the direct detection of coherent emission. This re- 
lies on the fact that an absorption coefficient is related 
to the concentration of absorbing centers in a simple and 
precise way. 

After discussing the principle of adiabatic spin refocus- 
ing in Section II, we present the absorbance detection of 
spin refocusing in Section III and IV. The experiment, 
described in Section V, is conducted in Tm 3+ :YAG, an 
actively investi gat ed system in the prospect of quantum 
storage 
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When spins interact with a spatially non-uniform static 
vertical magnetic field, the energy level spacing is inho- 
mogencously broadened. Then, excitation by a RF hor- 
izontal pulse gives rise to an oscillating magnetization 
that rapidly vanishes as the various spins precess at dif- 
ferent rates around the static field. It is well known that 
a 7r-pulse can rephase the spins and restore the oscil- 
lating horizontal magnetization. In order to uniformly 
excite all the spins, this pulse must be shorter than the 
inverse inhomogeneous broadening. This condition may 
entail inaccessible or unrealistic RF power requirements. 
Indeed, at given pulse area, the needed RF power varies 
as the square of the inverse pulse duration. By succes- 
sively driving the spins with different energy level spac- 
ings, frequency-chirped RF pulses may offer a less power- 
demanding way of excitation. This is precisely the way 
ARP works. 

Since the early days of Nuclear Magnetic Resonance 
(NMR) , ARP has been extensively used to flip vertically 
oriented spins, or equivalently, to perform a total transfer 
of population between two energy levels [2I( . This works 
as well on optical transitions [22[. Adiabatic passage is 
generally performed by sweeping the RF field frequency, 
while the static field remains constant. The spin is driven 
by an effective magnetic field B e / f whose horizontal and 
vertical components are repectively equal to the field am- 
plitude and proportional to the detuning. As the field fre- 
quency is swept, 'Beff rotates, flipping for instance from 
downward to upward orientation if the detuning is varied 
from negative to positive values. The spin adiabatically 
follows B e // provided the rotation speed is much smaller 
than the precession rate of the spin around B e y f . Vertical 
flip or population transfer, the most popular ARP appli- 
cation, does not preserve any phase-shift since the initial 
state is orthogonal to the final one. On the contrary, spin 
refocusing aims at making the final superposition state 
coincide with the initial one. The phase control is es- 
sential. A single ARP is not enough but a pair of ARP 
pulses can refocus the spins in any direction (23rf2o] |. 

The principle of spin refocusing by a pair of ARP pulses 
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FIG. 1: (color online) Time evolution of spin angular posi- 
tions in the horizontal plane. The red and blue lines corre- 
spond to two spins with different detunings. Initially aligned, 
the two spins return to their initially orientation at time T, 
after undergoing two ARP excitations centered at times T/4 
and 3T/4. The black dotted lines represent the locus of the 
flipping positions. 



is illustrated in Fig. [TJ The figure represents the evolu- 
tion of the phase, i.e. the angular position, of two spins in 
the horizontal plane. The two spins differ by their detun- 
ing with respect to the center loq of the inhomogeneously 
broadened energy spacing distribution. They are aligned 
at initial time t = 0. Then rotating at different rates 
around the vertical axis, they depart from each other. 
One applies a first ARP pulse whose frequency is swept 
at rate r through ujq at time T/4. The pulse gets resonant 
with a spin of frequency to at time T /A+(w— Wo)/r, when 
the phase equals rt(t — T/4). At that moment the spin 
vertical component is flipped and the horizontal angular 
position is reverted. In the figure, nipping is assumed 
to last much less than T/4. Then the spin phase con- 
tinues to evolve at the same rate until interaction with 
a second ARP pulse, identical to the first one, and cen- 
tered at time 3T/4. The interaction time is shifted by 
the same amount (w — ojo)/r with respect to 3T/4. This 
second pulse reverts again the phase and restores the ini- 
tial vertical spin component. Finally the spin are phased 
together again at time T. The locus of the flipping po- 
sitions in the phase/time representation is comprised of 
two quadratic segments defined by <j) = rt{t — T/4) and 
cf> = r(T -t)(t-3T/4). 

A simulation of spin rcfocusing by a pair of succes- 
sive ARPs is presented in Fig. [21 Each ARP lasts for 
T/2 =0.1ms, during which the RF field Rabi frequency is 
maintained fixed at 0.28MHz. The RF is swept at a con- 
stant rate of 0.04MHz//^s. The 2000 spins are distributed 
over a 1 MHz-wide gaussian inhomogencous distribution. 
At time t = 0, they are all directed along axis Ox. The 
figure illustrates the strong scattering of spin directions 
at times T/4, T/2 and 3T/4, strikingly contrasting with 
rcfocusing at time T. 




FIG. 2: (color online) Simulation of spin refocusing. The 
spins, distributed over a lMHz-wide inhomogeneous distribu- 
tion, are all directed along Ox at time t — 0. Their position 
on the unit sphere is represented at times (a), T/4 (b), T/2 
(c), 3T/4 (d) and T (e), where T = 0.2ms. Two ARPs are ap- 
plied successively, for 0.1ms each. They are frequency-chirped 
at the constant 0.04MHz//^-rate. They are swept through 
the center of the inhomogeneous distribution at times T/4 
and 3T/4 respectively. Their Rabi frequency is kept fixed at 
0.28MHz. 



III. MEASURING THE SPIN REFOCUSING 
EFFICIENCY 

Aiming at measuring the spin rcfocusing efficiency, we 
have to compare the final oriented spin concentration 
with the initial one. This information cannot be pro- 
vided by a standard spin echo experiment. Indeed, it is 
difficult to relate the echo intensity to the initial spin co- 
herence concentration. Usually an echo experiment only 
gives access to the signal intensity dependence on various 
parameters such as the pulse delay, temperature, exter- 
nal fields etc... Optical absorption is a simple tool to 
reach the desired quantity. 

Let the two states |1) and |2) of a / = 1/2 nuclear 
spin be connected by optical transitions to a common |e) 
excited state (see Fig (3]). The spins are first prepared 
in state |1) by optical pumping from state |2). In other 
words, all the spins are aligned vertically, with the same 
orientation. The transmitted optical intensity is origi- 
nally measured as I\ on transition |2) <— \e). As the 
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spins arc pumped from \2) to |1), the sample absorption 
vanishes and the transmitted intensity raises to the input 
value Iq. 
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FIG. 3: Optical probing of magnetic resonance. The upper 
electronic state |e) is optically connected to 2), a nuclear spin 
sublevel of the electronic ground state. As the |1) — |2) spin 
transition is excited by the RF field at Rabi frequency Q, the 
state 1 2) population is optically probed by absorption on the 
|2)-|e) transition. 

At time t = 0, one applies a RF tt/2 hard pulse to 
rotate the spins into the horizontal plane. Following the 
pulse, the vertical spin component drops to 0, as can be 
monitored through the absorption increase on transition 
|2) <— |e). States |1) and |2) are equally populated, which 
actually makes the absorption coefficient recover its origi- 
nal thermal equilibrium value. The transmitted intensity 
drops to I±. 

A tt/2 hard pulse, applied immediately after the first 
one but with opposite sign, should be able to rotate back 
the spins to their original orientation. Instead, one lets 
the spins process freely after the first tt/2 pulse. Ini- 
tially aligned, the spins begin to depart from each other 
in the horizontal plane, until they are evenly oriented in 
all directions. This occurs in a time interval equal to 
the inverse inhomogeneous width. After the spin direc- 
tion is randomized, excitation by a tt/2 hard pulse is no 
longer able to rotate back the spins to state |1) and to 
make the absorption vanish on transition |2) <— |e). The 
transmitted intensity remains equal to I\. 

To make the absorption sensitive again to a tt/2 pulse, 
one refocuses the spins with a pair of ARP full pulses 
that shine the sample at times T/4 and 3T/4. At time 
T, the spins are aligned again. A tt/2 pulse is used to 
rotate back the spins to vertical position. However, the 
horizontal spin component may have relaxed during the 
time interval T. If the spins do not totally return to 
state |1), absorption does not vanish. With exponential 
relaxation, at rate 7, the final absorption on transition 
1 2) <s— |e) is reduced by factor 1 — e _7T with respect to 
the original value, in the medium at thermal equilibrium. 



IV. REPLACING tt/2-HARD PULSES BY 
ADIABATIC HALF PASSAGE 



Resorting to tt/2 hard pulses is somehow inconsistent 
with previous assumptions. As already discussed, the 



uniform excitation of the inhomogeneously broadened 
spin transition would demand very brief hard pulses with 
unacceptably high power. Adiabatic passage may offer 
an alternative way to rotate the spin into the horizontal 
plane. Let us consider the transformation operated by 
an adiabatic half passage (AHP) process [2J], an ARP 
that is interrupted at time t = 0, when the pulse fre- 
quency is scanned through ujq. We restrict the discussion 
to the most simple case where the RF field amplitude 
is kept constant over the pulse duration, although more 
sophisticated adiabatic schemes, such as the complex hy- 
perbolic secant, turn out to be more popular in NMR 
and optical (2(| [27} applications. 

According to Appendix IA 21 the AHP turns the ini- 
tially vertical spin M(A, —00) = (0, 0, 1) into: 



M(A,0) 



1 



V^ 2 + A 2 



(—ft, 0, A) 



(1) 



where the Rabi frequency fi = 7b Ti is expressed in terms 
of the RF field constant amplitude B and the gyromag- 
netic factor jb- The spins with large negative A detuning 
undergo an upside down flip, while those with large pos- 
itive detuning preserve their initial orientation. Hence 
M Z (A, 0) is an odd function of A. Spins with |A| < fl 
are rotated into the horizontal plane, irrespective of the 
detuning sign. 

Let the spins be distributed over a Ao-broad inhomo- 
geneous width, according to an even distribution G(A), 
with unit normalization /G(A)dA = 1. Then the ab- 
sorption coefficient on transition |2) <— |e) reads as 



(f) = ao / G(A)[l-M,(A,t)]dA 



(2) 



At equilibrium, M Z (A, —00) = 0, and a(— 00) = a^. Af- 
ter pumping from |2) to |1), a(t) vanishes. The AHP 
makes a(t) change back to a$ at t = 0, exactly as a tt/2 
hard pulse would do, the transmitted intensity returning 
to h. 

However, the vertical spin component only vanishes 
on average. When the refocusing AFP pulses excite the 
medium, the spins with different A values do not flip si- 
multaneously, which is reflected by a transient variation 
of absorption. Let the spins evolve freely until the first 
AFP pulse, centered at time T/4, with A T/4 >> 1. As 
far as absorption is concerned, inhomogeneous dephas- 
ing enables us to ignore the conversion of the horizontal 
components into the vertical one. Indeed this contribu- 
tion vanishes when averaged over A. According to Ap- 
pendix [X3l M z (A,t) varies as: 
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The resulting variation of a(t), given by Eq. ([2]), carries 
information on G(A). After reaching a minimum: 
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at T/4, a(t) returns to ao, while the transmitted inten- 
sity reaches a maximum J 2 before dropping back to I\ . In 
a similar way, the transmitted intensity decreases during 
the second AFP pulse at 3T/4. 

At time T the spins are refocused as: 



M(A, T) = 



1 



Vn 2 + A 2 



(-tte^ T ,0,A) 



(5) 



where the transverse relaxation has been accounted for. 
An AHP is used again to rotate the spins back to vertical 
orientation. The pulse frequency is swept with opposite 
rate with respect to the first AHP, starting at frequency 
Wo at time T. The final M z component reads as: 



M(A,+oo) 



n 2 c -< T + a 2 



(6) 



dA 



(7) 



n 2 + a 2 ' 

which leads to the final absorption: 

The final intensity If(T) is related to io and I2 by: 

ln(J / (T)/I 2 ) = e^ T ln(/ // 2 ) (8) 

Should horizontal refocusing be imperfect, e _7T should 
be replaced by rye _7T where rj < 1. Hence, the refocusing 
efficiency can be defined as: 



V = e^ T bx(I f (T)/I 2 )/]n(I /I 2 ) 



(9) 



Provided 77 does not vary with T, Eqs. (|SJ) and ([9]) can 
be recombined into: 

ln(I f (T)/I f (oo)) =e-^ T - T °hn(I f (T )/I f (cx>)), (10) 

Hence, the variations of In (I/(T)) with T give access to 
the 7 decay rate, irrespective of 77 value. 



V. EXPERIMENTAL 

The experiments are carried out in a 0,1 % at. 
Tm 3+ :YAG crystal. An external magnetic field lifts the 
nuclear spin degeneracy. The resulting 4-level structure 
is comprised of two ground states and two excited states. 
We choose the same field orientation as in Ref. [JJI . 

The optical aspects of the setup have been described 
extensively in Refs.Jlfl EH- Basically, the light beam, 
emerging from an extended cavity diode laser, is am- 
plitude and phase-shaped by acousto-optics modulators, 
driven by a high sample-rate arbitrary waveform gener- 
ator (AWG). The crystal is cooled-down to 1,7 K in a 
liquid Helium cryostat. The static magnetic field is gen- 
erated by superconductive coils and is set to about 0,5 
T, which leads to a ground state splitting of 15,7 MHz. 
The spin transition is resonantly driven by a RF mag- 
netic field. This excitation is conveyed to the crystal by 



a 10-turn, 20 mm long, 10mm in diameter, coil oriented 
along the light pulse wave vector. The crystal sits at the 
coil center. The RF signal, generated by the AWG, is fed 
to the coil through a 500 W amplifier. 

First we calibrate the Rabi frequency of the spin tran- 
sition, as a function of the voltage applied to the coil. 
To this end, we optically detect the Rabi oscillations of 
the spin states. Let's go back to the system described in 
FiglU After optically pumping the ions in state |1), we 
apply the RF excitation and detect the state oscillations 
by monitoring the optical absorption on transition |2)- 
|e) with the help of a weak probe beam. Provided the 
probe intensity is low enough, the absorption coefficient 
is proportional to the population of |2) and oscillates at 
the Rabi frequency as observed in Fig. SJa). They are 
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FIG. 4: (color online) Nutation experiment, (a) time evolu- 
tion of the normalized transmission through the crystal. The 
transmitted optical intensity, equal to L when both ground 
state sublevels are equally populated, raises to Io when all 
spins are prepared in state When the RF magnetic-field 
is applied, the transmitted intensity oscillates (red line). The 
RF field is turned on at t=0. The peak-to-peak voltage at 
the coil is set to 120 V. (b) Rabi frequency as a function of 
the applied peak-to-peak voltage. 

damped because the spins are not uniformly coupled to 
the RF field. This partly results from the inhomoge- 
neous broadening of the spin transition, and from the A 
detuning dependence of the V^ 2 + A 2 effective Rabi fre- 
quency. This also reflects the non-uniformity of the RF 
field over the crystal, given the finite length of the coil 
and the field distorsions by metallic parts of the cryostat. 
The excitation non-uniformity also explains the limited 
amplitude of the oscillations in Fig. Eta). Total inver- 
sion is not observed. From the data in Fig. HJb), O is 
measured to vary linearly with the peak-to-peak voltage, 
measured at the coil input, at the rate of 2.4±0.1 kHz/V. 

Next we turn to the adiabatic spin refocusing experi- 
ment. We use the most simple adiabatic passage proce- 
dure where a fixed-amplitude field is frequency-chirped 
at constant rate r. As pointed out in Section [Hi the 
precession speed of the spin around the driving vector, 
Q, shall be much larger than the rotation speed of the 
driving vector, r/ft. In addition, the duration of the 
pulse shall by far exceed the time Ao/r needed to scan 
the inhomogencous width Ao. Finally, since an ARP is 
a coherent process, the flipping time Q/r shall be much 
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FIG. 5: (color online) Measurement of In (//(T)/7/(oo)) as 
a function of delay T. Exponential decay fit leads to 7" 1 = 
0.33 ± 0.05ms. The dashed (red) line corresponds to exponen- 
tial decay at rate 7 = 3.0ms _1 . The shaded area is inacces- 
sible. Indeed, the pair of refocusing AFPs prevents us from 
reducing T under 200/iS. However, refocusing is not needed 
if one applies the second AHP just after the first one. From 
the measurement of In (I f (T) / 1 f (00)) at T = (red dot), one 
infers that If(T) decays at rate smaller than 7 = 3.0ms -1 
from T = to T = 0.2ms. 



smaller than the spin coherence lifetime. The adiabatic 
pulse (AFP) features are summarized in Tab. |TJ They 
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TABLE I: Features of the adiabatic pulses used for spin refo- 
cusing 

are expected to satisfy the three adiabatic passage con- 
ditions. The spin coherence lifetime can be deduced 
from the T-dependcnce of the probe beam transmitted 
intensity If(T), measured at the end of the second AHP 
(see Eq. (fTO)) ). According to the experimental data, dis- 
played in Fig. [5j In (If (T)/If (00)) decays exponentially 
as expected, which leads to a spin coherence lifetime of 
7 -1 = 0.33 ± 0.05 ms. This result is consistent with a 
previous Raman echo measurement (l7j . The exponen- 
tial decay is also consistent with photon echo experiments 
that were performed in a Tm 3+ : YAG crystal at the same 
concentration, and in the same temperature conditions 
where the relaxation of both optical and spin coherences 
is dominated by spin-spin interaction p8j . However the 
spin coherence appears to decay at slower rate between 
T = and T = 0.2 ms. This effect deserves further 
investigation. 

The refocusing efficiency can now be measured. The 
time-evolution of the probe beam transmitted intensity is 
shown in Fig.[6l together with the RF adiabatic sequence. 
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In accordance with the storyline of Section IIV1 the in- 
coming probe intensity Iq is totally transmitted at the 



FIG. 6: (color online) spin refocusing experiment in a 0.1 at. 
% Tm 3+ :YAG crystal. Upper box: amplitude and frequency 
variations of the RF driving field. Lower box: transmitted op- 
tical intensity. After optical transfer to state |2) the medium is 
transparent (lo). The first AHP rotates the spins and restore 
the equilibrium transmission (Ii). The transmitted intensity 
raises to a maximum I2 during the AFP and behaves in the 
opposite way during the second AFP. At the end of the second 
AHP the transmitted intensity reaches its final If(T) value. 



beginning, through an intially transparent material. The 
first AHP excitation rotates the spins into the horizontal 
plane, which makes the transmitted intensity drop to I\ 
at time 0. Then, during the AFP refocusing steps, the 
transmitted intensity remains close to I\ , except around 
times T/4 and 3T/4. It first exhibits a bump, reaching 
the maximum value I2 at times T/4, then an anti-bump 
at time 3T/4. Those variations correspond to the flip of 
the vertical spin component. Indeed only the spins at 
wo have been perfectly rotated into the horizontal plane 
by the AHP. Off- resonant spins, distributed over the |1)- 
|2) transition inhomogencous width Ao, are only partially 
rotated and keep a vertical component. A magnified view 
of the bump at T/4 is shown in Fig. [7] We determine Ao 
(full width at half maximum) by fitting the bump profile 
with Eq. ©. We obtain A = 0.5 ± 0.1 MHz, which is 
consistent with previous measurements [lj} . 

After interaction with the two AFPs, the spins are ex- 
pected to be aligned together at time T. At that moment, 
after changing the chirp rate sign, we use a second AHP 
to rotate back the spins to their initial vertical orien- 
tation. Perfect refocusing, with no decoherence, should 
make the transmitted intensity return to Iq. The mea- 
sured transmitted intensity is denoted If(T). Substitut- 
ing the Fig. [5] experimental data and 7 -1 = 0.33ms into 
Eq. ©; we obtain 77 = 1.6 ± 0.3. This unphysical value, 
much larger than the unit upper boundary, suggests that 
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FIG. 7: (color online) transmission bump at T/4. The trans- 
mitted optical intensity temporal profile, as derived from 
Eq. Q, is fitted to the experimental data, assuming a gaus- 
sian inhomogeneous distribution. The RF frequency is swept 
at a rate of 40MHz/ms, and the Rabi frequency has been 
measured to be Q/(2ir) =0.28MHz. The dotted (blue), 
dashed(red) and dotted-dashed(magenta) lines respectively 
correspond to an inhomogheneous bandwidth ( full width at 
half maximum) of 0.3, 0.5 and 0.7 MHz. 



the spin coherence relaxation might be slower than ex- 
pected. We suspect that the 0.2 ms interval between the 
adiabatic passages is not much larger than the reservoir 
fluctuation correlation time. As a consequence, the spin 
decoherence might be partially suppressed by the two 
successive AFPs [2!| [3(J. This is confirmed by experi- 
mental data in Fig. [5] 



VI. CONCLUSION 

We have used optical detection to investigate the 
spin refocusing capabilities of adiabatic pulses in a 
Tm 3+ :YAG crystal. As in any conventional spin echo 
experiment we have measured the transverse relaxation 
time of the Tm 3+ ion spin coherence in the electronic 
ground state. Moreover we have shown that the opti- 
cal absorbance procedure gives direct access to the spin 
coherence concentration. As a consequence, one can 
measure the spin refocusing efficiency. Spin decoherence 
seems to be slowed down on time scales shorter than, 
or on the order of, 0.2 ms. If confirmed, this system 
might offer an interesting test bed for extending the Carr- 
Purcell Meiboom-Gill [3, HJ method to adiabatic pas- 
sages. 

This work is supported by the European Commission 
through the FP7 QuRep project, and by the national 
grant ANR-09-BLAN-0333-03 



Appendix A: Adiabatic Rapid Passage description 
1. ARP operator 

This very standard calculation is presented for the sake 
of completeness. In addition we want to point out the 



importance of a common reference frame in a problem 
where the medium undergoes successive excitations by 
RF fields Bi(t). Those fields are swept through cj at 
different times tj. They can be expressed as: 



Bi(t) = B cos [cvot + 4>{t — ti)} 



(AT) 



Instead of changing immediately to the frame rotating at 
the instantaneous frequency of the field, we first make a 
change to the common reference frame rotating at fre- 
quency ujq. In this frame, the rotating wave approxima- 
tion of the Bloch equation reads as: 



d f M x 



d t M y 



-AM y - fl sin [4>{t - ti)] M z 
AM X + fl cos [<j>(t - U)] M z 
d t M z = n sin [4>{t - U)] M x - cos [<j>(t - U)] M y 

(A.2) 

where A = CJ12 — ujq and the Rabi frequency Q = jbB 
is expressed in terms of the RF field constant amplitude 
B and the gyromagnetic factor jb- One performs two 
successive changes of frame. First one transforms the 
spin coordinates according to the rotation M = M' 
where: 



cos[(j>(t - U)] - sin [<f>{t - U)] 
sin [<f>(t- ti)] cos [<j>{t-ti)] 
1, 



(A.3) 



The new frame rotates in concert with the frequency 
swept field, and the Bloch equations are changed into: 



d f M x , = - 



A-(j>{t-ti) 



d t M v , = 
d t M z , = 



A — cj)(t- U) 



My, 

M x > + OAL ( A -4) 



-QM y i 



To align the z" axis along the driving vector, one per- 
forms a second transformation that is defined by the ro- 
tation M 1 = R { 2 l) M", where: 



4 i] (t) 



cosG^ -sinOW 

1 
sinOW cos9W 



(A.5) 



where: 

sin6 (i) = Sl/toeff, 
and 



cos 



A-4>{t-u) /n eff , 



n eff (t) 



o 2 



A-<p(t-U) 



In this frame the Bloch equation reads as: 



d t M x „-M z „Q^ 
d t M y n = <>, f/ .\/,.. 

atM z ,, + M x „ew 



-n eff My» 



(A.6) 



7 



Assuming adiabatic conditions are satisfied, one can ne- 
glect Ma»eW and M X „B^ , which leads to: 



d t M x 



-n ef fM y >, 



d t M„„ = n eff M x „ 



d t M z n = 



(A.7) 



The solution of this equation can be expressed as the 
following propagator: 



' cos$W -sin$W 6 s 
— sin<I>W cos$W o 
V 1, 



(A.8) 



where: $^ = f t , fl e ff(t")dt". Finally, in the common 
frame rotating at u>o, the effect of an ARP on the spin 
coordinates can be described by the propagator: 



Rf ] (t)R { 2 > {t)U {l) (t <- 1 7 ) R [ i> (t')R$> (0 (A.9) 



AHP effect 



The effect of the first AHP is described by (t <- 
-co). One easily verifies that R 2 °\—oo) = 11. As a 
consequence, the operator 



t/ (0) (io 



go) 4 0) (-oo)#T(-oo) 



(0)/ 



has no effect on the initially vertical spin M(A, — oo) 
(0,0, 1). One can also check that i?i (to) = H~ Hence 



M(A,t ) = 4 0) (*o)M(A,-oo) 



where, according to Eq. (|A.5|) . 



-^2 0) (*o) — 



eff 



a o -fr 
o 1 o 
n o a 



(A.10) 



(A.11) 



and O e // = Vfi 2 + A 2 . The action of the second AHP 
can be analyzed in a similar way. 



comes into coincidence with cjq. Since the spins have 
precessed for a time much longer than the inverse in- 
homogeneous width, their A-averaged horizontal compo- 
nent vanishes. Therefore, the vertical component evolves 
alone and, quite in the same way as in Appendix IA 21 
the AFP effect is completely described by R 2 (t). More 
specifically, 



M z (A,t) = M z (A,t ) 



A-^t-tO 



n 2 + A-0(t-ti) 



(A.12) 



where M z (A,i ) is given by Eg. lA.lOl 



Appendix B: Spin refocusing 

The RF fields £>i(t) and B^(t), defined according to 
Eq. (|A.1|) . are centered at times ti = t + T/4 and 
t;s = to + 3T/4. Interacting with those fields, a spin 
precessing in a static field undergoes two successive iden- 
tical ARPs. The A detuning is << ujq. According to 
Eq. (|A.9[) . propagation from to to to +T can be expressed 
as^ (3) (t 4 <- t 2 )V (1) {t 2 <- to), where t 2 = t +T/2. Since 
the two ARPs are identical, t/ (3) (t 4 <- t 2 ) = C/ (1) (t 2 <- 
to), Rf\t ) = Rf\t 2 ), and Rf\t 2 ) = Rf\t 4 ), accord- 
ing to the above definitions (see Eqs (|A.3I) and (|A.8[0 . 
With the additional assumption <p(t) = 4>(—t), all the rel- 



evant R^ (tfc) matrices are identical. It is also assumed 
that il << |A| + |0(t o - ti)| so that each ARP can be 
completed. Then R 2 (to) = 1 and: 



R ( 2 \t2) 



provided <p(t - ti) < 0. If <j>(t - h) > 0, R^\t ) and 
i?^(t2) expressions should be swapped. Combining all 
these properties, one easily shows that 




3. AFP effect on the vertical spin component 



^ (3) (t 4 ^t 2 )F (1) (t2^to) = lL 



We consider the vertical spin component evolution in 
the vicinity of time ti, when the first AFP frequency 



At to + T, the three spin components return to their to 
initial values. 
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